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Abstract. Spatial relation is one of the basic topics of geographic information 
systems (GISs), and it has been widely applied in data organization and spatial 
analysis. It seemed whether single topological, metric or directional relation can 
hardly represent complex spatial relations in some cases. Different topological 
relations may result in the same expression by some topological relation models, 
among which 9-lntersection Model (91 M) and its extended models (e.g. DE+9IM) are 
well known. However, we are used to describe spatial relation in several parts to 
make it more accurate, e.g., "Go north from here and along the lake about 200 
meters, and then you will arrive." I n this paper, the detailed and integrated model of 
the topological, metric and directional relation is proposed, to describe the complex 
spatial relations between region and region, line and line, line and region more 
exactly. We use one dimensional matrix to represent topological and directional 
relations, and two dimensional matrixes to represent metric relation. It's better to 
answer the question, e.g., "How the road crosses the park?" We also give an 
application finally, to calculate the actual area of the land use parcels, from which the 
area of related linear features is to be subtracted, according to their topological 
relations in details. 
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1. Introduction 

Spatial relation is one of the basic topics of geographic information systems (GISs), 
and it has been widely applied in data organization and spatial analysis. Spatial 
relation includes topological relation, metric relation and directional relation. Many 
scholars have focused on the progresses of spatial relation researches, among which 
the topological relations are studied mostly. Egenhofer et al. (1991, 1993) described 
9-lntersection model (91 M), and other extended models by 91 M such as DE+9IM, 
CBM (Clementini 1993, 1995), V9IM (Chen et al. 2000) , Three-Valued 
9-lntersection Model (Kurata 2009a, 2010). Randell et al. (1992) proposed Region 
Connection Calculus (RCC), which uses connection C(x, y) to describe relations 
between objects. Chang et al. (1987) took 2D- String to detect topological relations 
between objects according to their projections on X and Y axis. About metric 
relations, Egenhofer et al. (1998), Shariff et al. (1998) used splitting to determine 
how a region's interior, boundary, and exterior are divided by a line's interior, 
boundary, and vice versa; closeness measures to describe how far apart disjoint parts 
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are. About directional relation, the main methods as follows: Minimum Boundary 
Rectangle (M BR) (Papadias et al. 1996), CONE (Haar 1976), 2D- String (Chang et al. 
1987) and etc. 

And in recent years, Duboisset et al. (2007), Egenhofer(2009), Schneider et al. 
(2004,2006), Praingetal. (2008), Mckenneyet al. (2006,2007),Kurata(2008),Zhilin 
et al. (2006), Kimfung et al. (2007) discussed the complex topological relations 
between composite features. 

There also are some scholars who integrated topological, metric and directional 
relation (Sharma 1996, Hong 1996, GODOY et al. 2004, Cicerone et al. 2003, 
Schwering 2007 ,Dubeet al. 2009, Kurata 2009b), but they mostly paid attention to 
interactive reasoning of these three relations instead of their integrated 
representation. Actually, in practice, we often describe these three relations together, 
e.g., "How long and where does this tunnel cross the hill? How long does the road 
pass along the lake?" We are used to describe spatial relations in several parts to 
make it more accurate, but not simply or equivocally. We may say, eg., "Go north 
from here and along the lake about 200 meters, then you will arrive." Therefore, it 
tends to be more in accordance with our lives when the topological relation, metric 
relation and directional relation are integrated together. 

Spatial relations between point and point, point and line, point and region are 
seemed to be simple, but those between region and region, line and line, line and 
region are more complex. So we here only aim at these complex relations, to make 
out their detailed and integrated representation, suppose our objects are polylines 
without arcs, but regions can be convex or concave. 

The remainder of this paper is structured as follows: Section 2 proposes one 
general integrated method to describe topological, metric and directional relation. 
Section 3 takes some examples of spatial relations between region/ region, line/ line, 
and line/ region to prove the method above. Section 4 shows the applications of our 
research, for an example that how to calculate the actual area of land use parcels 
according to the detailed topological relations. Section 5 draws some conclusions and 
discusses future works. 



2. Integrated representation of spatial relations 



If R (A, B) represents the composite spatial relation between A and B, it can be split 
into n parts of single spatial relations, which areR (Al Bi), R (A 2 , B 2 )...R (A n , B n ): 

"R(Ai,Bi)~ 



R (A, B) = 



R(A 2 ,B 2 ) 



(1) 



R(An, Bn) 

If we use T, M, Dto describe respectively the topological relation, metric relation, 
and the directional relation between A and B, then R (A, B) can be represented as 
follows: 
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R(A,B)=T(A,B) v M (A, B) v D(A,B) (2) 
Whereas T(Ai, Bi), T (A 2 , B 2 ),...,T(A n , B n ) are serial topological relations between 

(AM, (A 2 ,B 2 ) (An, B n ), and D(A L Bi), D(A 2 ,B 2 ) D(A n , B n ) are serial directional 

relations; suppose there are k kinds of metric relations(e.g., length, area) between A 
and B, recorded by Ml, M2,...,Mk, then we get Mh, Mk,.., Mln to describe the first 
metric relations between (Al Bi),(A 2 , B 2 ),..,(A n , B n ), and Mki, Mk 2 ,..., Mk n to describe 
the k kind metric relations, see Figure 1 
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Figure L Topological, metric, directional relations integrated model 

We use one dimensional matrixto describe topological and directional relations, 
and two dimensional matrixes to describe metric relation, as Eq. 3: 





"T(Ai,Bi) _ 




Mi(Ai,Bi) 


M 2 (Ai,Bi) • 


• Mk(Ai,Bi)" 




^D(Ai,Bi)" 


R (A, B) = 


T(A 2 ,B 2 ) 


V 


Mi(A2,B 2 ) 


M 2 (A 2 ,B 2 ) • 


• Mk(A 2 ,B 2 ) 


V 


D(A 2 ,B 2 ) 




T(An,Bn) 




Ml(An,Bn) 


M 2 (An,Bn) • 


• Mk(An, Bn) 




D(An, Bn) 



3. Examples 



3.1. Relations between region and region 

According to studies of Egenhofer, there are 8 basic topological relations between 
region and region, i.e., disjoint, meet, equal, inside, contains, covers, covered by, 
overlap. But in Figure2, the topological relation between A and B can hardly be 
described by"overlap". However, it can be split into two parts, to getR (A, B) =R (Al 
Bi) + R (A 2 , B 2 ). If only the intersection area is considered as their metric measure, 
then we get: 



R (A, B) = 



T(Ai,B0 
T(A2,B2) 



V 



Area(AinBi) 
Area(A2 n B2) 



v 



D(Ai,B0 
D(A2,B 2 ) 



(4) 
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I n which: T (A L Bi) =T (A 2 ,B2) = 



Suppose Area (A 2 nB 2 ) =100m 2 , D (A L Bi) =D (A 2 , B 2 ) =North. 
The natural language of spatial relation between A and B as: A is in the north of B; 
A is partly overlapped and adjoined by B, the intersection area of which is 100m 2 . 




Figure 2. Spatial relationship between complex regions 
3.2. Relations between line and line 

As shown in Figure3, it's really complex for the relation between A and B, because 
both intersection and overlay exist. Now this complex relation can be structured in 
three parts, i.e., R (A, B) = R (A L Bi) + R (A 2 , B 2 ) + R (A 3 , B 3 ). Length is the metric 
measure to illuminate how long these two lines overlap. Directional relation is not 
considered here, si nee people seldom care for it. 



R (A, B) = 



~T(Ai,Bi)~ 




Len(Ai n Bi) 




T(A 2 ,B2) 


V 


Len(A2 n B2) 


(5) 


T(A 3 ,B3) 




Len(A3 n B3) 





I n which: T (A L Bi) =T (A 2 ,B 2 )= 



$ $ -.0 
$ -.0 


T(A 2 ,B 2 ) = 




-10 -10 -10 




-^(j> —iif) -^(j> 



Suppose :Len(A L Bi)=lDOm, Len(A 2 ,B 2 ) =Len(A 3 ,B 3 )=0 

The natural language of spatial relation between A and B as: Line A is partly 
overlapped by B for 100m, and then disjointed from B, finally they intersect at one 
poi nt. 



Figure 3. Spatial relationship between complex polylines 



3.3. Relations between line and region 




(a) (b) (c) (d) 

Figure 4. Four different topological relationships between line and region but similarly 
represented by 91 M and DE+91 M 



As shown in Figure4, there are four different relations between a line and a region. It 
seems to lines all pass through the outside, boundary, and inner of regions, i.e., lines 
all "cross" regions; even more, the expressions of 91 M and DE+91 M are the same 
(Eq.6,Eq.7). 



T (A, L) = 



A°rU° A°C\dL A°f)L 
dA[\L° dAf}dL dA[M~ 
A~[~]L° A~ftdL A~C\L 



-\<j> —\<j> 
-\<j> —\<j> 



(6) 



Dim (A, L) = 



dim(A° fU°) 


d\m(A°[\dL) 


dim(A°f]L ~) 




1 - 2 


dim(dAf)L°) 


d\m{d A[\dL) 


dim(dAf)L~) 




1 1 


dim(A~ fU°) 


dim(^~ (1 dL) 


dim(A~f]L~) 




1 2 



(7) 



Indeed, there are obvious differences in the order, number, and length of the 
crossing sections. 

For anypolylineA and region B, when they intersect we get n nodes. So L may 
splitinto n-1 straight lines, recorded byl_i,L2,.., L n -i.Thus, there are only three classes 
for theinner lineLi : a)l_i° <= A" , ( i.e., the line holy outside of the region); b) l_i° <= d 
A , ( i.e., the line holy at the boundary of the region); c) l_i° cA°,( i.e., thelinein 
the region). Now it would not happen that the line partly I ay outside of the region and 
partly in the region. We can get 6 classes of topological relations, their semantics are: 
inside, meet, disjoin, out touch, in touch, overlap. 

Taking Figure4 (a) as an example, if L is a road, A is a park, then we get three 
parts after intersecting, suppose their lengths are 200m, 300m, 100m. Thus the 
topological relation between the road and park can be expressed by: R (A, L) = R (A, 
Li)+R(A, L 2 )+R(A, L 3 ). 



(8) 



The natural language of spatial relation between A and B as: This road firstly lies 
in the boundary of the park for 200m; then partly in the park for 300m; finally 
outside of the park for 100m. As a result, now we can know clearly "how does this 
road pass through the park (M ark and Egenhofer, 1994)?" 
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4. Applications 

In Chinese land use database, a continual region which has the same land use 
class is mapped as a parcel, and a zonal land use class is surveyed by its sketch line 
and mapped as a linear feature, such as a small road or channel within the farm. Let 
GS be the parcel's geometrical area, and LS be the total area of the linear feature 
within the parcel , then the parcel's actual area (AS) should be calculated as follows: 
AS=GS-LS. As shown in Figure5 (a), a farm is passed by a long channel, which has 
just a short width (e.g., 2 meter). But on the small scale map (such as 110000) the 
channel is represented as a single line I which has no area of geometry (Figure 5(b)). 
In the case, rather than the parcel's geometrical area, the actual area of parcel p 
should be reduced by the related area of channel within the farm, seetheshadow part 
in Figure5 (a). 



Figure 5. (a)A channel passed a farm i n reality; (b) the channel drawn asa linear feature on 
the small scale map (such as 1 10000) 

We distinguish the different topological and metric relations between the parcel 
and linear feature in details, so as to accurately calculate how much area of thelinear 
feature with in the parcel to be subtracted from (i.e. LS). 

a) If the part of linear feature is outside of the parcel, its area need not to be 
subtracted; 

b) If the part of linear feature is entirely within the parcel, its all area should be 
subtracted; 

c) I f the part of I i near feature is al ong the boundary of the parcel , its area shoul d 
be subtracted half. 

As shown in Figure5(b), thelinel can be divided by some parts I X 12 and etc. 12 is 
entirely within (but in touch) the parcel p, so all the area of 12 (120m 2 =60mx2m) has 
to be subtracted from p; 1 1 passes along (meet) the boundary of p, so half of the area 
of II (40m 2 =40mxlm) to be subtracted; the other part of I is outside of p, so needn't 
to be subtracted. 



We developed a program to calculate the parcels' actual area, which passed by 
some linear features, the attribute lists show the results of actual area, as shown in 
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Figure6. Figure7 shows the records of metric relations, according the detailed 
topological relations between the parcels and linear features, whereas DKH means 
the parcel ID, XZBH means the linear ID related to the parcel, XZCD and XZKD 
mean the length and width of the line, while DLMJ is the area of the line to be 
subtracted. 
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Figure 6. A program developed to calculate the actual area of land use parcels 
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Figure 7. The records of metric relations and detailed topological relations between the 
parcels and I i near features 



5. Summary 

Spatial relation isso much important for GISsand its applications in reality that 
more people and academic conferences pay attention to it. In earlier years, most 
studi es focused on si mpl e and perfect spati al features, but fel I short of compl ex 
objects and can hardly represent complex spatial relations. Someones attempted to 
propose more compl ex models but it made things even complicated. I n this paper, we 
suggest a serial model which iscomposed by some single expression (eg. 9IM) in 
details, and integrates topological, metric, directional relations together to describe 
more exactly. H owever, this model is unshaped at present and needs to be promoted 
in future. 
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